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ABSTRACT 


The rRNA gene of the Willowsia bimaculata (Entomobryidae) were determined and com- 
pared the D3 fragment of the 28S rDNA with those of other collembolan species for their 
phylogenetic relationships. Sequences were analyzed using maximum parsimony, maximum 
likelihood, and neighbor-joining methods. The inferred phylogenies support monophyly of the 
order Collembola. The present result agreed with the phylogeny of suborder Arthropleona 
divided into two superfamilies (Poduroidea and Entomobryoidea) as mentioned in conventional 
systematics and also when compared with analysis of data of their morphological and allozyme 
characters performed and published elsewhere by the present author. However, concerning the 
suborder Symphypleona phylogenetic position, MP and NJ analyses suggest possibility that the 
Symphypleona may be a sister of the superfamily Poduroidea. An overall review of the results 
shows that the fast evolving part of 28S rDNA are informative for determining affinities 
among close relatives such as the family level and partly genus level in Collembola. 


Key words: Willowsia bimaculata, Collembola, Іпзесіз, 285 rRNA gene, variable region, 
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INTRODUCTION 


Collembola are the most widespread terrestrial 
arthropods with largest range of ecological diversity. 
They are usually found in damp situations. Many 
species inhabit soil, leaf litter, decaying wood, fungi, 
mosses, caves, and even the surface of water. With the 
exception of a few species, Collembola are considered 
to be useful organismS whose role in organic decom- 
position is critical to the long-term well being of soils 
(Greenslade, 1991; Hopkin, 1997). 

The classification and phylogenctic relationships 
among groups of the order Collembola have been a 
subject of disagreement among various authors without 
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teaching a satisfactory consensus (Paclt, 1956; Salmon, 
1964; Moen and Ellis, 1984; Soto-Adames, 1996; Bel- 
linger et al., 2001). Comparative external morphology, 
however, is still the main approach used in collembolan 
taxonomy. The Collembola were traditionally divided 
into two suborders each characterized by a distinct 
body form, Arthropleona with elongate body and 
Symphypleona with globular body. The suborder 
Arthropleona was classified in two superfamilies, Podu- 
roidea (Poduromorpha) with prothorax free and Entomo- 
bryoidea (Entotomobryomorpha) with prothorax more or 
less reduced. 

Paclt (1956), however, was the first to consider Podu- 
ridae іп Poduromorpha and Actaletidae in Entomobryo- 
morpha as sister taxa based on the hypognathous 
position of the hcad and place these taxa as the sister 
group of Symphypleona which also share this character. 
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The sister group relationship of Poduridae and Actaletidae 
to Symphypleona was subsequently supported by 
Salmon (1964). He defined suborder Metaxypleona for 
Podura and Actaletes. Moen and Ellis (1984) have 
suggested that Arthropleona and Poduromorpha should 
be discarded from formal classification because their 
definition is based on plesiomorphic characters and they 
may not form natural group. They considered a sister 
group relationship between Poduridae and Symphypleona- 
Neelipleona. Bellinger ct al.(2001) divided the 
Collembola into five suborders, based on a linear 
phylogeny of the progressively increasing function of 
the furca. In the tree showing them, based on phylo- 
‘genies proposed by Cassagnau (1971), Massoud (1976), 
Moen and Ellis (1984), Bretfeld (1986), Fjellberg (1994) 
and Soto-Adames (1996), Poduromorpha branched first, 
Metaxypleona next and Neelipleona as a sister group of 
Entomobryomorpha and Symphypleona. 

During the past decade, the ficld of molecular 
systematics which used various molecular data to infer 
the phylogenetic relationships among taxa has been 
rapidly developed. Changes in gene sequences are 
thought to be largely unaffected by the environment. 
And the similarity of DNA sequences may accurately 
reflect relationships among taxa (Avise, 1994; Нап, 
2000; Jung, 2000). However a given gene may be 
informative at one taxonomic level, but uninformative 
at another (Hopkin, 1997). 

Ribosomal RNA (tRNA) genes evolve at a relatively 
slow rate, which can provide optimal levels of sequence 
variation among higher taxonomic categorics such as 
phyla, classes, and orders (Hwang ct al., 1998; Min et 
al., 1998). Recently, the fast cvolving parts of 18S and 
28S rDNA have been used for the phylogenetic studies 
among close relatives at lower categories such as the 
family level in tiger beetles (Volger et al., 1997) and 
the genus level in Collembola (Carapelli ct al., 1995). 

Unfortunately, the molecular approach is in an early 
stage in Collembola systematics. I previously deter- 
mined. the 18S rDNA sequences from some Collem- 
bolan species (Lee et al., 1995а; 1995b). I found that 
the 18S rDNA sequence analyses of Collembola give 
contradictory results when the two different methods 
are used. They did not support taxonomic decision at 
family level. Recently, Frati and Dell’Ampio (2000) 
used the nuclear 28S rDNA and the mitochondrial 
СОП gene to establish phylogenetic relationships 
among species of the family Neanuridae. Тһе D3-D5 
fragment of the tRNA gene was, however, very 
conserved and supplied little information of relation- 
ships at lower category. The phylogenetic reconstruction 


Park KH 


based on СОП gene was partly in accordance with 
morphological data. 

The sequence analysis of the ТЕМА genc have been 
attempted little in the species of the family Entomo- 
bryidae in contrast with the some information available 
on these genes in some other families. Therefore, the 
TDNA sequence of the Willowsia bimaculata in Entomo- 
bryidae was compared with the D3 fragment of 28S 
ТОМА of 15 collembolan species in five families 
reported in the GenBank. It was hoped that the highly 
variable region of 28S rDNA might provide usefulness 
for reconstructing phylogenetic relationships within and 
among five collembolan families. This is the first 
attempt that phylogentic relationships among the 
families in Collembola are inferred by comparing the 
fragment of the 28S rDNA. 


MATERIALS AND METHODS 


Genomic DNA extraction 

The preparation of the collembolan genomic DNA 
was carried out with slight modification of the method 
of Edwards et al. (1991). Fifty individuals of Collembola 
were crushed in а 1.5 ml microtube with liquid 
nitrogen. After adding 500 1 of DNA extraction buffer 
(100 mM Tris-HCI, рН 8.0; 80 mM EDTA, pH 8.0; 
160 mM Sucrose) to each tube, they were mixed by 
inverting the tubes several times. SDS (10%) and 
proteinase К (final concentration 15 mg/ml) were 
added, mixed, and incubated overnight at 37°C. Then 
equal volume of chloroform and phenol were added. 
The tube was spun at 3,000 rpm for 5 min and the 
supernatant was transferred to a new tube to be mixed 
with two volume of cold absolute ethanol to precipitate 
the genomic DNA by centrifugation for 10 min at 3,000 
rpm. The genomic DNA was dissolved with 20 1 of 
double-distilled water. 


Cloning and sequencing 

The sequences were obtained from cloned PCR 
fragments. A fragment of the rDNA cistron, including 
the ITS2, was amplified by polymerase chain reaction 
(PCR) and the primers were as follows: c58 (5’-GAA- 
CTGCAGGACACATGAACA-3’) апа с280 (5’-CTGA- 
CCAGGCATAGTTCACC-3’). I designed forward and 
reverse primer scquences in the conserved coding 
region of the 5.8S and 28S rRNA genes based on 
known sequences of insects. The PCR reaction mixture 
comprised of 50 ng templete DNA, 25 pmol of each 
primer, 4 zl of 2.5 mM each dNTP, 5 xl of 10x 
reaction buffer, 0.5 unit of Tag DNA polymerase 
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Table 1. List of Taxa Used for Phylogenetic Comparisons. 
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Таха Species Abbreviation Accesion Number 
Tnsecta 
Collemboia 
Arthropleona 
Poduroidea 
Poduridae Podura aquatica Poduaqua AF005468 
Neanuridae Friesenae Friesea grisea Friegri AJ251739 
Pseudachorutinac Anurida maritima Anurimar А1251738 
Neanurinae Bilobella aurantiaca Biloaura AJ251730 
Bilobella massoudi Bilomass AJ251732 
Neanura muscorum Neanmus AJ251733 
Entomobryoidca 
Isotomidae Tetracanthella sp. Tetracan AJ009851 
Isotomurus nebulosus Isomuneb AJ009850 
Isotomurus indipendente Іѕотиіпа Х84953 
Isotomurus ghibellinus Ізотив йі X84957 
Isotomurus palustris Іѕотираі! X84952 
Isotoma klovstadi Isotmakl Y11704 
Entomobryidae Orchesellinae Orchesella ranzii Orchranz А1009852 
Entomobryinae Willowsia bimaculata Willbima AY066017 
Symphypleona 
Sminthuridae Allacma fusca Allacfus X84956 
Sminthurus viridis Sminvir AJ251741 
Thysanura 
Lepismatidae Thermobia domestica Therdome AF370808 


(TaKaRa Ex Taq.). The PCR reaction was performed in 
a GeneAmp PCR System 2400 thermocycler (Perkin 
Elmer Co.). Thirty amplification cycles were performed, 
each involving denaturation at 94°C for 1 min, 
annealing at 60°C for 30 sec, and extention at 72°C for 
| min. The reaction was linked to a final polymeri- 
zation step of 7 min at 72°C. The amplificd products 
were extracted again with GeneClean П (Bio 101) and 
diluted in 20 yl of distilled water. The PCR products 
were cloned into the pGEM-T easy vector (Promega 
Со.) and transformed into Е coli JM109 cells. Double- 
stranded recombinant plasmid DNA was purified using 
a plasmid purification kit (Promega Co.). The amplified 
rDNA fragments sequenced with an ABI Prism 377 
automatic DNA sequencer using BigDye terminator 
cycle scquencing kit V2.0 (Perkin Elmer Co.) in Korea 
Basic Science Institute, Kwangju Branch. Both 
directions were fully sequenced. 


Data analysis 

The taxa used in this study are listed in Table 1. The 
sequence data of D3 stem of 28S rDNA in Willowsia 
bimaculata has been submitted to the GenBank 
(accession number: AY066017). The D3 fragments of 


the 28S rDNA nucleotide sequences were retrieved 
from the GenBank detabase with the aid of a BLAST 
network server (Altchul et al., 1997). The sequence of 
Thermobia domestica (Thysanura) was employed as an 
outgroup for all analyses. Nucleotide sequences of 17 
species were initially aligned using the Clustal W 
(Thomson et al., 1994) alignment program before 
performing the phylogenetic analysis. Sites with 
missing data or gaps were excluded from all analyses. 
After aligning, a set of sequences were confirmed to 
the format of the PHYLIP package ver. 3.5 (Felsenstein, 
1993). From the realigned sequences, 1000 bootstrap 
samples were generated using the SEQBOOT program. 
The resulting matrices were subsequently subjected to 
the DNAPARS, DNAML and NEIGHBOR programs 
employing the maximum parsimony, maximum likeli- 
hood and neighbor joining analyses. Distance matrices 
were calculated using the DNADIST program and 
nucleotide substitutions (excluding gaps) were estimated 
using the Kimura’s two parameter method (Kimura, 
1980). A maximum likelihood tree was inferred using a 
transition/transversion ratio of 2:0. The trees were 
reduced to a single consensus tree using the CON- 
SENSE program. Finally, a most probable phylogenctic 
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ATTGGGTGTTTACGAAACCTAAAGGCGTAATGAAAGTAAAGIGTTTTCG-ATTACTTCGGTGGTCGAGAACCGAGGGAAGATGGAGTTTG 89 
ATTGAGAATGTT--AAACTCAAAGGCGTAATGAAAGTGAAGA---TCCA-ATTTTATGGAT—TTGAG----GAGAGATCATATCTTTIG 78 
ATTGGAAATTACT-AAATCCAAAGGCATAATGAAAGTAAAGA-—-TCCA-ATTTTATGGAT--ATGAG-—-GAGAGATCATGGCTTTAG 79 
ATTGGGAATTATT-AAACCCAAAGGCATAA TGAAAGTAANGA——TCCA-ATTTTATGGAT—-ATGAG--—-GAGAGATCATGGCTTTAS 79 
ATTGGGAATTACT-AAACCCAAAGGCATAATGAAAGTAAAGA——-TCCA-ATTTTATGGAT--ATGAG----GAGAGATCATGGCTTTAA 79 
ATTGGGAAATACT-AAACCCAAAGGCATAATGAAAGTAAAGA-—-TCCA-ATTTTATGGAT--ATGAG-———-GAGAGATCATGGCTTTAA 79 
ATTGGGAATTATT-AAACCCATAGGCATAA TGAAAGTASAGA—-TCCA-ATTTTATGGAT—-TIGAG-—--GAGAGATCATGGCTITGG 79 
ATTGGAAATGT-—-AAACCCAAAGGCATAATGAAAGTAAAGAA—TTCG-ATTGC—GAAT--TIGAG-—-GAGAGATCAAGGGCAGAG 76 
ATTGGGAATGT—-AAACCCAAAGGCATAATGAAAGTAAAGAA—-TTCG-ATTGC—-GAAT--TTGAG--—-GAGAGATCAAGGGCAGAG 76 
ATTGGGAATGT---AACCCAAAGGCATAATGAAAGTAAAGAA--TTCG-ATTGC--GAAT--TTGAG-——GAGAGATCAAGGGCAGAG 76 
ATTGGGAATGT---AAACCCAAAGGCATAATGAAAGTAAAGAA-—T'TCG-ATTGC--GAGI—TTGAG--——GAGAGATCAATGATAGAG 76 
ATTGGGAATGT---AAACCCAAAGGCGGAATGAAAGTAAAGAA--TTCG-ATTGC—GAGT—TTGAG-—--GAGAGATCAATGATATAG 76 
ATTGGGAAAGT-—-AAACCCAAAGGCATAATGAAAGTAAAGAA——TICG-ATTGC—-GAAT—TTGAG----GAGAGATCATGGAT-GAG 75 
ATTGGGAATGA-—- AAACCCAAAGGCGTAATGAAAGTGAAGA-—-TTCG-ATTGC—GiAAT—TTGAG----GAGAGATACCGTAGCGGT 75 
ATTGGGAATGA--AAACCCAAAGGCGTAATGAAAGTGAAGA——TTCG-ATTGC—-GAAT—TTGAG----GAGAGATCCTGGCICTGY 75 
ATTGAGAATGTA—-AAACTCAAAGGCATAA TGAAAGTGAAGAT—TCGT-CTTGTACGGAT-—-TGAG--—GAGAGATCATGTCAATAC 78 
ATTGAGAATGTA—-AAACTCAAAGGCATAATGAAAGTGAAGA-—TCCGTCTTGCACGGAT-—-TGAG-—GAGAGATCATGGCAATAC 78 


180 
TCAAGITIGTTCTGATCCITTTCGGAGGTGAGGTCGGGCTTGGCAACTTGATCCGAGTCAATTTGGCAAAGGATACTCCCGCACTCCCAG 179 
e ‘SAGTA--TTCCGTA----AGGT--—--TT----ATTAAAAAAA--—-GTATG--------------CGCATCTCCTG 124 
ATCT-——-----TAATA--TTCTGAA----AAG'---TT--—AT TAAGTAAA---—G TATG----------—--CGCATCTCCTG 127 
GICT------—-TGATA—-TTCCGTA----AGGT-—---T --—-ATTGAGTAAA---—GTATG----———--CGCATCTCCTG 127 
GTCT-----—-TGATA--TTCCGTA----AGGT--—--TT-—-AT TGAGTAAA-----GTATG-——----------CGCATCTCCIG 127 
GICT------—-TGATA--TTCCGTA——-AGasT--—-TT —ATTGAGTAAA-----GTATG----------—-CGCATCTCCTG 127 
A-CT ТОАТА--ТАТІСАЛ----АбАТ----- ТТ----АТТСАОТЛАА-----СТА10----------СОСАТСТССТ6 126 
А--------00006--ТТТСОТА----АбАС-----ТС---АСТОТСТСТ0---ТТТІ6--------<ССАТОТССІТО 121 
һ---------10066--ТТІСОТА----АСАС----ТС---АСТСІСІСІб----ТІТІ6----------СЧОСАТСТССТО 121 
А ---16006--ТТТСОТА----АСАС-----ТС----АСТӨТСІСІ0-----ТТТТ0------------СССАТСІССТО 121 
А-----------1(М45%--ТІССОТА----АСбТ-----ІС----ААТСІСІСІС----ТТТТЕ —————-CGCATCICCIG 121 
A ———TGGGG—TTCTGTA——AAGA-----TC--—-ACTGTCTATG-----TTATG-—----------—CGCATCTCCTG 121 
Mga ——-TGGGA--TTCCGCA-—-AGTT----- TT ——-CCGCTCT-1G-----TTATG--—-----------CGCATCTCCTG 120 
GTGG-----—--TGGC-—-TTCCGTA---AGGT----- T--—-GTT'GCCACAGC----TT1GGT-—-—------—-CGCATCTCCG; 122 
= —-—-TGGCA--TTCCGTA—-——- A(T ----TT ——-GCTGGGCAGG-----GTGGG--------------CGCATCTCCGG 120 
T-—-----——- TGGCA—-TICOGTA-—-AG-T- ——- TT -—-GCTGGGSTAAATT—-GTATG-------—-—-CGCATCTCCIG 121 
T———--—GGGCA—-TTCCTAA-—--AGGT-----TT- GC TCGGTACTT I —-GTATG---_---- EACE G 125 


GGCGTCTCGTTCTCATTGCGAGATGAGGCGCACCTAGAGCGTACACGTTGGGACCCGAAAGATGGTGAACTATGCCCGGCCAG 252 
GACGTCTT-AAC-CTTCGA-AGGTGGGGCGTACCTAGAGCATACACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 204 
GATGTCTT-AAC-CTCIGA-AGGTAGGGCA TACCTAGAGCATACACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 207 
GATGTCTT-AAC-CTCTGA-AGGTAGGGCATACCTAGAGCATACACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 207 
GATGICTT-AAC-CTCTGA-AGGTAGGGCATACCTAGAGCATACACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 207 
GATGTCTT-AAC-CTCTGA-AGGTAGGGCATACCGAGAGCATACACGTIGGGACCCGAAAGATGGIGAACTATGCCTGGICAG 207 
САССТТТТ-АЛС-СТТТСА-АССТАСССССТАССТЛСАССАТАСАССТТСССАССССАААСАТССТСААСТАТСССТССТСАС 206 
GACGTCTA-TAC-CICTGA-AGGTTGGGCGTACCAAGAGCATATACGTIGGGACCCGAAAGATGGTGAACTATGCCTGGICAG 201 
GACGTCTA-TAC-CTCTGA-AGGTTGGGCGTACCAAGAGCATATACGT TGGGACCCGAAAGGTGGTGAACTATGCCTGATCAG 201 
GACGTCTA-TGC-CTCIGA-AGGTTGGGCGTACCAAGAGCATATACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 201 
GACGTCTA-TAC-CTCTGA-AGGT ['GGGCGTACCAAGAGCATATACGTTGGGACCCGAAAGATGGTGAACTATGCCTIGGTCAG 201 
GACGTCTA-TAC-CTCTGA-AGGTTGGGCGTACCAAGAGCATATACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 201 
GACGTCTAATAC-CTCTGA-AGGTTGGGCGTACCAAGAGCATACACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 201 
GACGTCTA-TAC-CTCTGA-AGGTAGGGCGTACCTAGAGCGTACACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 202 
GACGTCTA-TAC-CICTGA-AGGTGGGGCGTACCTAGAGCGTACACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGICAG 200 
GACGTCTA~TAC-CTCTGA-AGGTTGGGCGTACCTAGAGCATACACGTTGGGACCCGAAAGATGGTGAACTATGCCTGGTCAG 204 
САССТСТА-ТАС-СТСТСА-АССТОСССССТАССТАСАССАТАСАССТТСССЛССССАААСАТССТСААСТАТСССТОСТСАС 205 


Figure 1. Sequence alignment of the D3 fragment of 28S rDNAs from 16 species of Collembola including 
Thermobia domestica (Thysanura) as an outgroup (- indicates alignment gap). 


tree in each analysis was produced using the Treev32. bimaculata was determined and aligned with the D3 


fragment of 28S rDNA sequences of 15 Collembolan 


RESULTS AND DISCUSSION species in five families and one Thysanuran specics as 


an outgroup. The nucleotide sequences of the 16 


The rDNA nucleotide sequence of the Willowsia specics except Willowsia bimaculata were obtained 
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from the data of GenBank (Table 1). The total aligned 
set of nucleotides of the 17 species is shown in Figure 
1. The aligned sequence length over all taxa is 263 bp. 
The multiple alignment of these 17 taxa showed that 
the sequence of Thysanuran species was longer than 
those of Collembolan species. In overviewing 17 sequen- 
cies, 1 bp insertion at 50 (Т) site апа 47 bp deletions at 
14 (б), 69-72 (AACC), 95-102 (СТТТОТТС), 108-109 
(СТ), 117-120 (GGTG), 125-129 (CGGGC), 132-135 
(GGCA), 148-150 (TCA), 156-169 (GCAAAGGATAC- 
TCC), 193 (Т) апа 200 (G) sites were noted in those of 
Collembolan species. 

In order to avoid the potential risks and inconsi- 
stencies of a single method or data set, the phylo- 
genetic reconstructions were based on three different 
approaches, methods of maximum parsimony (MP), 
maximum likelihood (ML) and neighbor joining (NJ). 
The inferred phylogenies support monophyly of the 
Order Collembola (bootstrap value: 100%). In the 
phyletic trees of MP (Fig. 2) and ML (Fig. 3) the orders 
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of branch within each family were same, but the 
position of Symphypleona were inconsistent in the two 
methods. In the MP tree, Collembola were subdivided 
into two branches: Subtribe 1 clustered with Sminthu- 
ridae which includes Allacma fusca and Sminthurus 
viridis, Poduridae which includes Podura aquatica, and 
Neanuridae which includes Friesea grisea, Bilobella 
aurantiaca, Bilobella massoudi, Neanura muscorum and 
Anurida maritima, whcreas Subtribe 2 clustered with 
Isotomidae which includes /sotoma Klovstadi, Tetracan- 
thella sp., and 4 species of Zsotomurus (1 nebulosus, І. 
indipendent, I. ghibellinus and 1. palustris), and Ento- 
mobryidae which includes Orchesella ranzii and 
Willowsia bimaculata. In the ML tree, Collembola were 
also divided into two subtrees but the position of 
Symphypleona was different from that of MP tree: One 
contains Poduridae and Neanuridae, the other clade 
includes Isotomidae and Entomobryidae plus Symphy- 
pleona. The phylogenetic reconstruction based on NJ 
(Fig. 4) was partly in accordance with MP and ML 
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Figure 2. The most parsimonious tree using DNAPARS in PHYLIP package. The numbers 
above the branches are confidence levels of consensus tree from 1000 bootstrap replicates 
(A: Poduroidea, B: Entomobryoidea, A+B: Arthropleona, C: Symphypleona, A+B+C: 


Collembola). 
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Figure 3. The maximum likelihood trec using DNAML іп PHYLIP package. The numbers 
above the branches are confidence levels of consensus tree from 1000 bootstrap replicates 
(А: Poduroidea, В: Entomobryoidca, A+B: Arthropleona, С: Symphypleona, A+B+C: 


Collembola). 


analyses, but it was disaccordant for the placement of 
some species. Isotoma klovstadi was а sister taxon to 
the Entomobryoidea clade, instead of ingroup to the 
taxa of Isotomidae. However, because of the lower 
bootstrap value (43%), the position of it is not obvious. 

The clustering pattern of present results agreed with 
the phylogeny of suborder Arthropleona divided into 
two superfamilies of Poduroidea and Entomobryoidea 
as mentioned in conventional systematics and also 
when compared with analysis of data of their 
morphological and allozyme characters (Park et al., 
1996) performed and published elsewhere by the 
present auther. In cladistic analysis of morphological 
characters and allozyme analysis using four species 
representing each family of the superfamily Poduroidea 
and three species of families of the superfamily Ento- 
mobryoidea, Arthropleona were subdivided into two 
branches. However, the 185 rDNA sequence analysis 
yielded only a limited support to this scheme. The NJ 
tree based on 18S rDNA sequence data supported 


clearly Poduroidea grouping but not the Entomobry- 
oidea one. The MP method, however, did not support 
monophyly of the superfamily Poduroidea but suggested 
the single grouping of three families of Entomobryoidea. 
The 18S rDNA sequence analyses gave contradictory 
results when the two different methods were used. 
These results showed the primary sequence from 18S 
TRNA gene is not suitable for determining taxonomic 
status at lower categories such as the family level in 
Collembola, or suggest that two superfamilis of Collem- 
bola may not form natural groups. In present study, the 
relationship among two superfamily groups is obvious, 
and the clustering pattern of the present molecular dada 
analyses agreed with the traditional view of these two 
superfamilies. The taxonomic position of Poduridae is 
by no mean certain as it shares some characteristics 
with the Arthropleona and others with Symphypleona 
and Neelipleona (Hopkin, 1997). Frish (1978) found that 
the heart of Poduromorpha, including Podura, has six 
ostia, and in Symphypleona the ostia are reduced to 
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Figure 4. The neighbor joinning tree using NEIGHBOR in PHYLIP package. The numbers 
above the branches are confidence levels of consensus iree from 1000 bootstrap replicates 
(A: Poduroidea, B: Entomobryoidea, A+B: Arthropleona, C: Symphypleona, A+B+C: Col- 


lembola). 


three. According to Hemmer (1990), the presence of a 
haploid chromosomes number in P. aquatica of п = 11 
compared to typically n=5 in Symphypleona makes a 
close relationship between these two taxa questionable 
(Hopkin, 1997). Based on these Frish (1978) and 
Hemmer (1990) concluded that the characters invoked 
by Paelt (1956) and Salmon (1964) to relate Poduridae 
and Symphypleona represented convergence and not 
teal homology. The present results based on 28S rDNA 
analysis, сусп though other species of Podura were not 
included, strongly suggests that this species show greater 
affinity to Neanuridae (Poduromorpha, Arthropleona) 
than to Symphypleona. This result supports the pers- 
pective of Bellinger et al. (2001) based on a linear 
phylogeny of the progressively increasing function of 
the furca. 

Although the Collembola was traditionally classified 
in two suborders, Arthropleona and Symphypleona, the 
phylogeny have been a subject of disagreement among 
various authors. Soto-Adames (1996) divided it into 


four groups (Poduromorpha, Entomobryornorpha, Symphy- 
plona, and Neelipleona), while Bellinger et al. (2001) 
classified it in five groups (including Metaxypleona). 
And several major questions concern the relatedness 
among suborders have been present. Some authors re- 
garded the Symphypleona as more closely related with 
the Poduridae (Paclt, 1956; Salmon, 1964; Moen & 
Ellis, 1984), whereas others considered closer to the 
Entomobryomorpha (Soto-Adames, 1996; Bellinger et 
al., 2001). We obtained contradictory results in the 
position of the suborder Symphypleona when the three 
different methods are used. Close affinity of Symphy- 
pleona and Poduridae resulting from the MP and NJ 
methods partly support the systems of Paclt (1956), 
Salmon (1964) and Moen and Ellis (1984). The relation- 
ship of them is strongly supported by bootstrap values 
(100% and 84%). However, the affinity of Symphypleona 
and Entomobryidae resulting from the ML тпсіһоа 
supports the systems of Soto-Adames (1996) and Bellinger 
ct al. (2001). However, because of the lower bootstrap 
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value (47%) of Entomobryidae and Sminthuridae, the 
telationship among them is not obvious. Therefore, the 
position of the Symphypleona is unstable throughout 
the three methods examined and requires more 
extensive examination to confirm its exact taxonomic 
position. However, MP and NJ analyses suggest 
possibility that Symphypleona may be a sister of the 
Poduromorpha. 

These results show the molecular evidence of highly 
variable region of 285 rRNA gene sequence аге 
informative for determining affinities among the family 
level in Collembola. However, the relationships of 
subfamilies or genus within Neanuridae are not clear 
and the arrangement presented here may not be a 
natural one. Particularly the relationships among the 
Pseudachorutinae and Neanurinae were unceriain. In the 
Isotomidae clade, /sotoma branched first апа Tetracan- 
thella next as a sister group to the taxa of Jsotomurus. 
The. relationships of genus within Isotomidae аге 
obvious. These results show the molecular evidence of 
the fast evolving part of 28S rRNA gene sequence 
might partly provide usefulness for determining 
affinities among close relatives such as the genus level 
in Collembola. 

It is universally accepted that the class Collembola is 
monophyletic group. However, classification schemes 
still use Linnaean categories which may be difficult to 
adapt to a stictly cladistic analysis of evolutionary 
relationships. The classification should be regarded for 
the moment as a convenient arrangement rather than 
the most parsimonious solution to the problem of how 
to subdivide the order Collembola (Hopkin, 1997). 
Several attempts have been made to reconstruct the 
evolution of Collemobola (Maynard, 1951; Yosii, 1961; 
Cassagnau, 1971; Uchida, 1971; Massoud, 1976). 
Nevertheless, a convincing phylogenetic tree has not yet 
been established. Obviously I need more sequences data 
of other collembolan and related species for the full 
picture of the evolution of Collembola. 
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